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I. INTRODUCTION 

In the last twenty years, Chiral Perturbation Theory 
(ChPT) Jl]-||] has emerged as a powerful tool to describe 
the interactions of the lightest mesons. These particles 
are considerably lighter than the rest of the hadrons, 
which is nowadays understood as a consequence of the 
spontaneous breaking of the SU(3)l X SU(3)r chiral 
symmetry down to SU(3)l+r that would be present in 
QCD if the three lightest quarks were massless. In such 
case, the light mesons would correspond to the massless 
Goldstone bosons associated to the spontaneous chiral 
symmetry breaking. Of course, quarks are not mass- 
less, but their masses are so small compared to the typ- 
ical hadronic scales, 0(1 GeV), that their explicit sym- 
metry breaking effect also translates into a small mass 
for the lightest mesons, which become pseudo-Goldstone 
bosons. Hence, the three pions correspond to the pseudo- 
Goldstone bosons of the SU(2) spontaneous breaking 
that would occur if only the u and d quarks were massless, 
which is a remarkably good approximation. Similarly, the 
meson octet formed by the pions, the kaons and the eta 
can be identified with the eight pseudo-Goldstone bosons 
associated to the SU(3) breaking, when the s quark is 
also included. 

The low energy interactions of pions, kaons and the eta 
can be described in terms of an effective Lagrangian that 
follows the SU(3)l x SU(3)r ~ > SU(3)l+r spontaneous 
symmetry breaking pattern. If we do not include any 
additional field apart from the pseudo-Goldstone bosons, 
this description will only be valid for energies much be- 
low the scale where new states appear. That is, the ef- 
fective ChPT Lagrangian provides just a low energy de- 
scription. As a consequence we can organize all the pos- 
sible terms that respect the symmetry requirements in a 
derivative (and mass) expansion. Therefore, any ampli- 
tude is obtained as a perturbative expansion in powers 
of the external momenta and the quark masses. The 



importance of this formalism is that the theory is renor- 
malizable and predictive, in the following sense: all loop 
divergences appearing at a given order in the expansion 
can be absorbed by a finite number of counterterms, or 
low energy constants, that appear in the Lagrangian at 
that very same order. Thus, order by order, the theory 
is finite and depends on a few parameters that can be 
determined experimentally. Once these parameters are 
known, any other calculation at that order becomes a 
prediction. Basically, these are the main ideas underly- 
ing ChPT, which has proven very successful to describe 
the low energy hadron phenomenology (for reviews see 

§)• 

Despite the success of this approach, it is unfortunately 
limited to low energies (usually, less than 500 MeV). That 
is the reason why, over the last few years, there has been 
a growing interest in extending the applicability range of 
the chiral expansion to higher energies. Of course, this 
requires the use of non-perturbative methods to improve 
the high energy behavior of ChPT amplitudes. These 
methods include the explicit introduction of heavier res- 
onant states in the Lagrangian resummation of di- 
agrams in a Lippmann-Schwinger or Bethe-Salpeter ap- 
proach or other methods that unitarize the ampli- 
tudes like the Inverse Amplitude Method (IAM) f,§. 
The latter has been generalized to allow for a coupled 
channel formalism || , yielding a successful description of 
the meson- meson scattering amplitudes up to 1.2 GeV, 
and even generating dynamically seven light resonances. 

In principle, these methods recover at low energies 
the good properties of ChPT, since they use part of 
the perturbative information. However, it should be 
noted that, so far, the full results to one loop for all 
the meson-meson scattering processes were not available 
in the literature. At present only the irn — > nn flTo| , 
Kit — > Kir |u^, t/tt — > r/ir [jl(J and the two independent 
K+K- -> K+K-, K+K- -> K°K° glj amplitudes 
have been obtained in the SU(3) ChPT framework, al- 
though with different procedures and notations. As a 
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consequence, the IAM has only been applied rigorously 
to the 7T7T, KK final states, whereas for a complete treat- 
ment of the whole low energy meson-meson scattering 
additional approximations had to be done In partic- 
ular, the lowest order expansion could not be recovered 
complete up to 0(p 4 ) thus spoiling the scattering lengths 
and, in addition, it was not possible to compare directly 
with the low energy parameters of standard ChPT in di- 
mensional regularization and the MS — 1 scheme. 

In this work, we have calculated all the meson-meson 
scattering amplitudes at one loop in ChPT. There are 
three amplitudes that have never appeared published in 
the literature: Kr\ — > Kr], r\r\ — ► r/rj and Kit — > Kr\. The 
other five have been recalculated independently and all 
of them are given together in a unified notation, ensuring 
exact perturbative unitarity and also correcting previous 
misprints. Then, we have applied the coupled channel 
IAM to describe the whole meson-meson scattering be- 
low 1.2 GeV, including low energy data like scattering 
lengths. This new calculation allows for a direct com- 
parison with the standard low energy constants of ChPT 
and that is why we have made a considerable effort in es- 
timating the uncertainties in all our results, which are in 
very good agreement with the present determinations ob- 
tained from low energy data without unitarization. The 
main differences of this work with are that we con- 
sider the full one-loop results for the amplitudes, ensur- 
ing their finiteness and scale independence in dimensional 
regularization, we take into account the new processes 
mentioned above and we are able to describe more ac- 
curately the low energy region. This had already been 
achieved for the irir KK system only in |ll[ , but here we 
complete this task for the whole meson-meson scattering. 

The paper is organized as follows. In section || we re- 
view the main features concerning the meson-meson scat- 
tering calculations at one-loop in ChPT. The final results 
for the amplitudes have been collected in Appendix B due 
to their length. The definition of partial waves and uni- 
tarity is discussed in section |l|, and the IAM is presented 
in section |y|. In section |v|, we review the availa ble d ata 



on meson-meson scattering. In sections VI and VII we 
first use the IAM with present determinations of the low 
energy constants and next we make a fit to the data com- 
mented in_section Our conclusions are summarized in 
section VIII. Apart from the amplitudes in Appendix B, 



we have also collected some useful formulae in Appendix 
A. 



II. MESON MESON SCATTERING AT ONE 
LOOP 

The lowest order Lagrangian for SU(3) Chiral Pertur- 
bation Theory is: 

C 2 = ^(d^d^U + M (U + U r )), (1) 



where /o is the pion decay constant in the SU(3) chi- 
ral limit and the angular brackets stand for the trace 
of the 3x3 matrices. The matrix U collects the 
pseudo-Goldstone boson fields tt, K, r\ through U (<E>) = 
exp(iv / 2 < I > //o), where 
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and Mq is the tree level mass matrix. Throughout this 
paper we will be assuming the isospin limit, so that Mq 
is given by 
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As a matter of fact, from these definitions, it can be 
easily seen that the tree level masses satisfy the Gell- 
Mann-Okubo relation [HJ: AM 2 K - M 2 „ - 3M 2 ^ = 0, 
that will be very useful for simplifying the amplitudes. 

From the Lagrangian in Eq.(|l|), one can obtain the 
0(p 2 ) amplitudes just by calculating the corresponding 
tree level Feynman diagrams. In order to obtain the 
0(p ) contributions, one has to consider loop diagrams, 
whose generic topology is given in Fig.Q, which will gen- 
erate UV divergences. If loop integrals are regularized 
with dimensional regularization, which preserves the chi- 
ral symmetry constraints, the divergences can be reab- 
sorbed in the chiral parameters Li of the fourth order 
Lagrangian: 

C 4 = L^dJJ^U) 2 + L 2 (d fl tfd v U)(d»tfd v U) (4) 
+L 3 (d fl U^d fl Ud„U ii d l/ U) +L 4 (d^d tl U)(U^Mo + M^U) 
+L 5 {d„U^d»U{U+M + M+U)) + L 6 {U^M + M+U) 2 
+L 7 (f/ t M - Mlu) 2 + L s {MlUMlU + U^MqU^Mq), 

where the terms which couple to external sources, like 
gauge fields, are omitted ^||. The Li constants are re- 
lated with the renormalized LK/i) generically as Li = 
LK/J,) + Ti\ [| where jj, is the renormalization scale, 



A 



d-4 r 



16tt 2 



l 



i(log47T- 7 



1) 



(5) 



7 is the Euler constant and the Ti coefficients can be 
found in ||. We remark that the L3 and L7 constants 
are not renormalized and are therefore scale independent, 
i.e, L 3 = T 7 = 0. 

Thus, up to fourth order one has to consider the tree 
level diagrams from 0(p 2 ) and 0(p 4 ), together with the 
one-loop diagrams in Fig.|l|. We stress that mass and 
wave function renormalizations should be accounted for 
to the same order. The latter are schematically repre- 
sented by the tadpole diagram (e) in Fig.[l| As custom- 
ary, we define the bare fields in terms of the renormalized 
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ones as 7r = Z^ 2 Tt ren and so on for the kaons and eta, 
so that scalar fields have finite canonical kinetic terms. 
Taking into account all the different contributions from 
diagrams of type (e) in Fig.[l] plus those tree level dia- 
grams coming from £4, one obtains: 

4 2 4A , n 9 \ 
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we will not need the mass renormalization of in any 
of our expressions. 

The meson decay constants are also modified to one 
loop, ft will be convenient for our purposes to write 
all the one-loop amplitudes in terms of a single decay 
constant, which we have chosen to be f n . For that reason 
and for an easier comparison with previous results in the 
literature, we also give here the result for the meson decay 
constants to one loop pj: 
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It is important to stress that both the physical masses 
in Eq. (||) and the decay constants in Eq. (|J) are finite and 
scale independent. 



with i = 7T, K, 77. 

Note that the wave function renormalization constants 
Zi contain a divergent part and they are scale dependent. 

As for the mass renormalizations, the physical pion 
and kaon masses are given in terms of the tree level ones 
as 1: 
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According to the chiral power counting, we have to 
use Eq.(^|) and Eq.(|J) only in the tree level part of the 
amplitudes. In fact, the mass renormalization Eq.(||) af- 
fects only the mass terms coming from the Lagrangian in 
Eq.(|l]) and not the masses coming from the kinematics 
of the corresponding process. As it will be seen below, 





d e 

FIG. 1. Generic One- loop Feynman diagrams that have to 
be evaluated in meson-meson scattering. 



Therefore, the one-loop ChPT scattering amplitude 
(renormalized and scale independent) for a given process 
will have the generic form: 

T(s, t, u) = T 2 {s, t, u) + T£°\s, t, u) + T 4 "™(s, t, u) (10) 

where s,t,u are the Mandelstam variables. Here, T2 is 
the tree level contribution from the Lagrangian in Eq. ([TJ) , 
whereas Tf° contains the fourth-order terms which are 
polynomials in s,t,u. Those polynomials have four pos- 
sible origins: tree level terms from the Lagrangian in 
Eq.(Q) proportional to L[, other polynomial terms pro- 
portional to Li with i — 4. ..8 coming from the mass and 
decay constant renormalization in eqs.(||) and (Q), terms 
proportional to /jj coming from tadpole diagrams ((d) 
and (e) in Figure [l]) and finally pure polynomial fourth- 
order terms which stem from our parametrization of the 
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one-loop functions (see Appendix A). Let us remark that, 
for technical reasons explained in section III.B, we have 
chosen to write all our amplitudes in terms of /„• only 
since, using eqs.(fpis), Jk and f v can be expressed in 
terms of f n , L\ and L\. In addition, T" ni stands for 
the contribution of diagrams (a) , (b) and (c) in Figure 
These contributions not only contain the imaginary parts 
required by unitarity but also yield the correct analytic 
structure for the perturbative amplitudes, as it will be 
discussed below. We remark that all the terms in T% nl 
will be proportional to the J and J functions defined in 
Appendix [A|. 

Using crossing symmetry it is not difficult to see that 
there are only eight independent meson-meson ampli- 
tudes. We have calculated these amplitudes to one loop 
in 577(3) ChPT. They are given in Appendix || Three of 
these amplitudes had not been calculated before, namely 
K°i] — > K°rj, r\r\ — > r\r\ and Kr\ — > K vr . For the rest, we 
have checked that our amplitudes coincide with previous 
results [lOillJ] up to differences in notation and different 
simplification schemes, equivalent up to 0(p 6 ). In par- 
ticular, since we are interested in the "exact" form of 
perturbative unitarity (see below), we have written our 
final results in terms of a single pion decay constant, f v , 
and we have used the Gell-Mann-Okubo relation taking 
care of preserving exact perturbative unitarity. Further- 
more, we have explicitly checked that all the amplitudes 
remain finite and scale independent. 

Finally, we wish to add a remark about 17 — 17' mixing, 
since the physical r\ is indeed a mixture of the U(3) octet 
and singlet pseudoscalars, whereas in this work we are 
only using the standard SU(3) ChPT. One may wonder 
then if our description of the r\ is just that of the pseu- 
doscalar octet component, since in this Lagrangian the 
singlet field is not an explicit degree of freedom. How- 
ever, it has been shown fLq ] that the standard framework 
results from an expansion in powers of the inverse powers 
of the "topological susceptibility" of the complete C/(3) 
Lagrangian. In that context the rj' is considered as a 
massive state (that is why it does not count as an explicit 
degree of freedom) but the singlet component generates 
a correction to L7. Note that, indeed, the mass of the r\ 
contains an L7 contribution, and that is why we can use 
M v in eq.(|J) with its physical value, whereas the Mq v is 
the one satisfying the Gell-Mann-Okubo relation exactly. 
Therefore, our approach can be understood as the low- 
est order approximation to the n — rf mixing problem, 
where all the effects of the mixing appear only through 
L7. Since we will only compare with data in states with 
one 77 at most, and below 1200 MeV, our results seem 
to suggest that this approximation, although somewhat 
crude, is enough with the present status of the experi- 
mental data. Indeed, we will see that the values that we 
obtain for L7 are in perfect agreement with those given 
in the literature (and this comparison can now be done 
because we have the complete one loop amplitudes renor- 
malized in the standard way). 



III. PARTIAL WAVES AND UNITARITY 

A. Partial waves 

Let us denote by T* b (s) the partial wave for the pro- 
cess a — ► b, i.e, the projection of the amplitude for that 
process with given total isospin I and angular momen- 
tum J. That is, if T^ b (s,t,u) is the isospin combination 
with total isospin 7, one has 
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dxPj(x)T* b (s,t(s,x),u(s,x)) (11) 



where t(s, x), u(s, x) are given by the kinematics of the 
process a — ► b with x — cos 9, the scattering angle in the 
center of mass frame. 

Note that we are normalizing the partial waves includ- 
ing a factor N, such that N — 2 if all the particles in the 
process are identical and N = 1 otherwise. Recall that, 
since we are working in the isospin limit, the three pions 
are considered as identical, so that N = 2 only for the 
7T7T — ► 7T7T and r\r\ — > T)T) processes. 

We shall comment now on the T* b amplitudes for every 
possible process involving 7r, K, rj. Using crossing symme- 
try and assuming isospin symmetry exactly, we will de- 
termine the number of independent amplitudes for each 
process. The discussion is general and there is no need 
to invoke ChPT, although we will refer to the results for 
the amplitudes in Appendix [b|, which are the one-loop 
ChPT results. 

• 7T7T — > 7T7T scattering: There is only one independent 
amplitude, so that one has 

T°(s, t, u) = 3 T(s, t, u) + T{t 1 s, u) + T(u, t, s), 
T 1 (s, t, u) = T{t 1 s, u) — T(u, t, s), 
T 2 (s, t, u) = T{t 1 s, u) + T(u, t, s), 

where T(s,t,u) is the 7r + 7r~ — > 7r°7r° amplitude. 
At one-loop in ChPT it is given in Appendix 
Eq.(||). 



• Kit — > Kn scattering: Crossing symmetry allows 
us to write the I = 1/2 in terms of the 7 = 3/2 one 

as 

T 1/2 (s, t, u) = i [3 T 3 / 2 (u, t, s) - T 3 / 2 (s, t, u) 

(12) 



Here, T 3 / 2 (s,i,u) is the K+ir+ 7T+7T+ ampli- 
tude, whose expres sion at one-loop within ChPT 
corresponds to Eq.(B5). 



• KK — > KK scattering: We can write the isospin 
amplitudes as 
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T°(s,t,u) = T ch (s,t,u) + T neu (s,t,u), 

T x (s, t, u) = T ch (s, t, u) - T„ eu (s, t, u), (13) 

where T c h and T neu are, respectively, the ampli- 
tudes for the processes K + K~ — > K + K~ and 
K°K° — > ivf+ii' . Their ex pres sions to one loop 
correspond to Eqs.(B7) and ( |B8|) , respectively. 



• _?T7A" — > ttt: scattering: In this case, one has 

T °( S) t,u) = ^- [T^ 2 (u, s, t) + T 3 / 2 (t, s, u) 
T l {s, *,u) = -j= [T 3/2 (u, s, t) - T 3 / 2 (i, a, u)] , (14) 

where T 3 / 2 (s,t,u) is the iY + 7r+ — > if + 7r+ ampli- 
tud e, g iven in Appendix [b| for one-loop ChPT, 
Eq.@. 



• Krj — > if 77 scattering: This is a pure 1=1/2 pro- 
cess. The one-loop amplitude can be read directly 
from Eq.( |B2| ) 

• KK — > 777/ scattering: This is an / = process that 
using crossing symmetry can be obtained from the 
previous amplitude as follows: 

• i^?7 — ► scattering: This is also an 1=1/2 
process, whose amplitude correctly normalized, is 

T l / 2 {s,t,u) = -V3T ROv ^ R0 ^ (s,t,u), (16) 

where the one-loop exp ression for K°r/ — > K°ir° 
can be found in Eq.(jB3]) 



• KK — > m] scattering: This is a I — 1 process re- 
lated to the K°r) — > K°ir a amplitude by crossing 
symmetry, i.e: 

^(SjtjU) = -V2T R o v ^ R07r o(t,s,u). (17) 



ttt] — > 7r?7 scattering: This is a pure 1=1 isospin 
amplitude whose one-l oop ChPT expression can be 
read directly from Eq.(B6). 



• 7T7T — > r/rj scattering Now 1 = and the amplitude 
is obtained from the previous one by crossing, as 

T„o„o_+ m (8, t, u) = T„o v _>„o v (t, s, u). (18) 



• 7777 — > 7777 scattering: Here, 1 = and the corre- 
sponding one-loop amplitude can also be read di- 
rectly from Eq.(Bl). 



In this paper we will be interested in the case when there 
are several coupled states for a given choice of 7, J, i.e: 
the coupled channel case. In particular, with the above 
normalization, the relationship between the T-matrix ele- 
ments T~i and the S'-matrix ones is given for two coupled 
channels (a, b = 1,2) by 



S n = 1 + 2io- 1 T n , 
S22 = 1 + 2i (J2T22, 



5*12 = S21 = 2i y/ (71(72 Tl2, 



(19) 
(20) 
(21) 



where the I J superscripts have been suppressed to ease 
the notation and we have used that due to time reversal 
invariance, 7y = Tji. Here, cr; = 2qi/^fs where is 
the center of mass momentum in the state i. Note that 
Oi is nothing but the phase space of that state at y/s. 
In the 1 = channel above the r/rj threshold we will 
use the corresponding generalization in the case of three 
channels. 



B. Unitarity 

The S matrix should be unitary, i.e, SS^ = 1. In case 
there is only one state available, that means that S can 
be parametrized in terms of a single observable, which is 
customarily chosen as the phase shift. For the case of two 
channels, the elements are organized in a unitary 2x2 
matrix, containing only three independent parameters. 
We will follow the standard parametrization: 



S = 



iy/l - ife l ^ +s ^ 



where the Si are the phase shifts and 77 is the inelasticity. 

The unitarity relation translates into relations for the 
elements of the T matrix of a particularly simple form 
for the partial waves. For instance, if there is only one 
possible state, " 1" , for a given choice of I, J, the partial 
wave T11 satisfies Eq.(|l9|), so that unitarity means 



ImTn = (7i I Tn 



ImTiI 1 



(23) 



In principle, the above equation only holds above 
threshold up to the energy where another state, "2", is 
physically accessible. If there are two states available, 
then the T matrix elements satisfy 

ImTn =(7i(s)|Tn| 2 + (72|Ti2| 2 , 
lm J 12 = (7i in J 12 + 02 I12 J 22 1 
ImT 22 = (7i I T12I 2 + cr 2 I T 22 | 2 . 

In matrix form they read, 

ImT = TST* =*> ImT- 1 = (24) 

with 
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Tn T12 

T12 T; 
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(25) 



which allows for an straightforward generalization to the 
case of n accessible states by using n x n matrices. 

One must bear in mind that the unitarity relations 
imply that the partial waves are bounded as the energy 
increases. For instance, in the one channel case, from 
Eq.(p3|) we can write: 



T - sinS P *s 
in — e 



(26) 



where S is the phase of tn. 

Note that all the unitarity relations, Eqs.(p3|) and (24), 
are linear on the left hand side and quadratic on the right. 
As a consequence, if one calculates the amplitudes per- 
turbatively as truncated series in powers of an expansion 
parameter, say T = T 2 +T4 + . . ., the unitarity equations 
will never be satisfied exactly. In particular, for ChPT 
that means that unitarity can only be satisfied perturba- 
tively, i.e.: 



ImT 2 = 0, 

Im T4 = T 2 E T2 , 



(27) 



where the latter is only satisfied exactly if one is care- 
ful to express T4 in terms of masses and decay constants 
consistently with the choice made for T 2 . That has not 
always been the case in the literature and that is one of 
the reasons why we have recalculated some processes: all 
our results satisfy exact perturbative unitarity. Other- 
wise there are additional 0(p 6 ) terms in Eq. (p7f) . As we 
will sec below, this will be relevant to obtain a simple for- 
mula for the unitarized amplitudes. Our choice has been 
to rewrite all the fx and f v contained in the amplitudes 
in terms of f^, L\ and L\ using the relations in eq.(||). 

The deviations from Eq.(|24|) are more severe at high 
energies, and in particular in the resonance region, since 
unitarityimplies that the partial waves are bounded, 
see Eq.(|2q), which cannot be satisfied by a polynomial. 
Generically, in the resonance region, the unitarity bounds 
are saturated. If a polynomial is adjusted to saturate uni- 
tarity in a given region, in general, it will break the uni- 
tarity bound right afterward. Another way of putting it 
is that resonances are associated to poles in the complex 
plane, that will never be reproduced with polynomials. 

For all these reasons, if we are interested in extend- 
ing the good properties of ChPT to higher energies, we 
have to modify the amplitudes, imposing unitarity and 
a functional form that allows for poles in the complex 
plane. This will be achieved with the Inverse Amplitude 
Method. 



IV. THE COUPLED CHANNEL INVERSE 
AMPLITUDE METHOD 

As it can be seen from the unitarity condition in 
Eq. (p4[) , the imaginary part of the Inverse Amplitude 
is known exactly above the corresponding thresholds, 
namely, ImT -1 = — E. Indeed, any amplitude satisfying 
the unitarity constraint should have the following form: 



T= (RcT- 1 -iE) 1 



(28) 



Consequently, we should only have to calculate the real 
part of T _1 . As a matter of fact, many unitarization 
methods are just different approximations to ReT" 1 (see 
H for details). The idea behind the Inverse Amplitude 
Method (IAM) is to use the formula right above, but 
approximating ReT -1 with ChPT. Since we have T ~ 
T2 + T4 + .... Then 



T- 1 ~T 2 - 1 (1-T 4 T 2 - 1 + ...), 
ReT" 1 ~ T 2 _1 (l - (ReT 4 )T 2 - 1 H 



(29) 
(30) 



so that multiplying Eq.(p8|) by T 2 T 2 1 on the left and 
T 2 _1 T2 on the right, we find 



T ~ T 2 (T 2 - ReT 4 - i T 2 E r 2 )- 1 T 2 . 



(31) 



At this point, if the amplitude satisfy "exact perturbative 
unitarity", namely Eq. (p7j) , we can simplify the above 
equation to obtain a simple expression 



T ~ T 2 (T 2 - T i )~ 1 T 2 



(32) 



This is the generalization of the IAM to coupled channels. 
Note that this formula only ensures exact unitarity if T4 
satisfies "exact perturbative unitarity" . 

The IAM was first applied to just one elastic channel 
and it was able to reproduce well the nir and ttK scat- 
tering phase shifts below the KK and K r\ thresholds, re- 
spectively. In addition it was able to generate the a (now 
called / (400 - 1200)) the p and the K* resonances [§. 
Furthermore, it was shown how the formula for the one 
channel IAM can be justified in terms of dispersion rela- 
tions [H], which allowed for the analytic continuation to 
the complex plane and the identification of the pole as- 
sociated to each resonance in the second Riemann sheet. 

In view of Eq.(p2), it may seem necessary to know 
the complete 0{p r ) ChPT calculation of each one of 
the T matrix elements. Nevertheless, one could use a 
further approximation and calculate only the s-channcl 
loops (Fig. la), which are the only responsible for the uni- 
tarity cut and are supposed to dominate in the resonant 
region. This was the approach followed in Q, having 
in mind that the complete ChPT calculations were not 
available at that time for any meson-meson scattering 
two-channel matrix. The results were remarkable, repro- 
ducing up to 1.2 GeV seven (/, J) meson- meson scatter- 
ing channels (17 amplitudes), and even generating seven 
resonances. However, the fact that the s-channel loops 
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were regularized with a cutoff, together with the omis- 
sion of crossed loops and tadpoles, made impossible to 
compare the chiral parameters with those of standard 
ChPT (still, they had the correct order of magnitude, 
as expected). Besides, the low energy ChPT predictions 
were recovered only partially. This motivated the au- 
thors in to calculate the full 0(p 4 ) K + K~ -> K+R- 
and K + K~ — > K°K° amplitudes]], which allowed for the 



low KK threshold, the IAM 7T7T scattering element will 
have an additional spurious contribution from the imag- 
inary part of the KK scattering left cut, which extends 
up to sfs = 4(M|- — M%). This is a well known and 



unitarization with Eq.(f32| of the {I, J)=(0, 0) and (1, 1) 
channels. This approach yielded again a good high en- 
ergy description but also reproduced simultaneously the 
low energy mr scattering lengths. All these results were 
obtained with Li parameters compatible with those of 
standard ChPT 0. 

As we have seen in the previous section, we have cal- 
culated the last three independent 0(p 4 ) meson-meson 
scattering amplitudes that were still missing. They are 
given in Appendix B in an unified notation with the other 
five that we have recalculated independently (correcting 
some minor misprints in the literature). Therefore, we 
are now ready to unitarize the complete meson-meson 



scattering by means of Eq.(32) 



However, at this point we have to recall that for a given 
energy Eq. (|32|) has only been justified for a matrix whose 
dimension is exactly the number of states accessible at 
that energy. The reason is that the unitarity relation, 
Eq. (p4[) , increases its dimensionality each time we cross 
a new threshold. Thus, for instance, in 7T7t scattering, 
one should use the one dimensional IAM up to the KK 
threshold, then the two dimensional IAM, etc.... although 
this procedure yields discontinuities on each threshold, 
instead of a single analytic function. Another possibil- 
ity H is to use the IAM with the highest possible di- 
mensionality of the /, J channel for all energies []. This 
second possibility yields an analytic (and hence continu- 
ous) function, but it may not satisfy unitarity exactly at 
all energies, namely, when the number of opened chan- 
nels is smaller than the dimensionality of the IAM for- 
mula. Following with the ttit — KK example, if we use 
the 2-dimensional IAM formula, we will have exact uni- 
tarity ensured above the KK threshold, but not below. 
In particular, if we still use the 2-dimensional IAM be- 



*An erratum for these amplitudes has appeared published 
when preparing this work. The previous results and conclu- 
sions in jnj are nevertheless correct, since the errata did not 
affect the numerical calculations. We thank J. A. Oiler for 
discussions and for letting us check that their corrected am- 
plitudes coincide with ours. 

'As a technical remark, let us note that in such case, the 
IAM has to be rederived in terms of the partial waves T a b 
divided by the CM momenta of the initial and final states, to 
ensure that this new amplitudes are real at lowest order. From 
there the derivation follows the same steps, and we recover the 
very same Eq. (^) by multiplying by the initial and final state 
momenta in the end. 



lasting problem in the literature 14-16 that affects also 
other unitarization methods, like the K-matrix |^5J . As 
a matter of fact, several years ago |l4| it was suggested 
that the physical solution would probably be an interpo- 
lation between the two just mentioned approaches. How- 
ever, in the context of ChPT and the IAM, and for the 
7T7T — KK channels, it was found that the violations 
of unitarity are, generically, of the order of a few per- 
cent only. We have confirmed this result but now for the 
whole meson-meson scattering sector. Even the thresh- 
old parameters can be accurately reproduced, since they 
are defined through the real part of the amplitudes which 
are almost not affected by the spurious part. The origin 
of this problem is that the IAM in Eq.(|3^) mixes the 
left cuts of all the channels involved when performing 
the inverse of the T2 — T4 matrix. Thus, it is not able 
to reproduce the left-cut singularities correctly I|[l7| , al- 
though numerically their contribution is negligible when 
all the observables are expressed in terms of real parts of 
the amplitude, and taking into account the present status 
of the data and the uncertainties in the Li. 

In this paper we have chosen to show the second ap- 
proach, since the one-dimensional IAM has been thor- 
oughly studied in ||. Very recently there have been 
proposed dispersive approaches |l6| to circumvent this 
problem in the 7T7t KK system, but they involve the cal- 
culation of left-cut integrals that are hard to estimate 
theoretically. It would be interesting to have them ex- 
tended and related to the ChPT formalism, but that is 
beyond the scope of this work. The fact that we use 
the higher dimensional IAM formalism, which contains 
spurious cuts, does not allow for a clean continuation 
to the complex plane. Nevertheless, since poles associ- 
ated to resonances have already been found in the one- 
dimensional case J^] and in other approximated coupled 
channel IAM approaches , we leave their description for 
a generalized IAM approach with better analytic prop- 
erties [Q. In this work we will concentrate on physical 
s values, and the compatibility of the unitarized descrip- 
tion of resonances and low energy data with existing de- 
terminations of the chiral parameters. Nevertheless we 
will also show that this can also be achieved with the 
first, discontinuous, approach. 



V. MESON-MESON SCATTERING DATA 

Let us then comment on the data available for each 
channel: 

Channel (I,J) = (1,1) For the energies considered 
here, the two states that may appear in this channel 
are 7T7t and KK. In Figs. 2a and 3a, we plot the data 
on the 7T7T scattering phase shift obtained from [Il9| and 
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pO| , which correspond to the squares and triangles re- 
spectively. Let us remark that the first set of data points 
tends to be between two and three standard deviations 
higher than the second when the phase shift is higher 
than 90 degrees, and the other way around for smaller 
values of the phase shift (note that error bars are smaller 
than the data symbols) . Thus the data sets are not quite 
consistent with one another, which could be fixed with 
the addition of a systematic error of the order of a few 
percent. 

This channel is completely dominated by the tttt state 
and there is almost no inelasticity due to KK produc- 
tion below 1200 MeV. The (1 - 77^ ) /4 points from the 
inelasticity analysis given in ]^l] are shown in the lowest 
part of Figs. 2. d and 3.d. 

Channel (I,J) = (0,0) For this channel we may have 
up to three states, namely tttt, KK and r\r\. In this 
case, there are three observables with several sets of data, 
which, as can be seen in Figs.2.b,c and d, are somewhat 
incompatible between themselves when only considering 
the errors quoted in the experiments. Again, they be- 
come compatible if we assume a systematic error of a 
few percent. For the tttt scattering phase shift (Sqo, see 
Fig.2.b), the experimental data shown come from: dif- 
ferent analysis of the CERN-Munich Collaboration [22] 
(open square), as well as from jl9) (solid square), | 2*3| 1 
(solid triangle) and (solid circle). Concerning the 
tttt — > KK phase shift, the data in Figs. 2. c and 3.c cor- 
responds to: (2^] (solid triangle) and ]^5| (solid square) 
and they are reasonably compatible, mainly due to the 
large errors in the first set. Finally, we are also show- 
ing in Figs. 2. d and 3.d the data for (1 — ?7qo)/4, since it 
is customary to represent in that way the values of the 
inelasticity 7700 ■ The experimental results are rather con- 
fusing here, mainly up to 1100 MeV, due to problems in 
the normalization. From the data shown in the figure, we 
have only fitted to those coming from: (25) (solid square), 
p6| (solid triangle), ]^7j (open square) and |2^] (open cir- 
cle). There is a disagreement in the normalization with 
the data of ]^9| up to a factor of 2 (see for discus- 
sion). We have not included the latter in the fit, mostly 
because in the analysis of |2|] they neglect the unitarity 
constraint, which in our approach is satisfied exactly at 
those energies. 

Channel (I,J) = (2,0) There is only the tttt state and 
so we only display in Fig.2.e and 3.e the 620 phase shifts 
again from the CERN-Munich Collaboration [31| (open 
square) and the CERN-Saclay Collaboration 32] (solid 
triangle). 

Low energy K14 decay data. This reaction is partic- 
ularly important since it yields very precise information 
on the Sqq — 5±i combination of tttt scattering phase shifts 
at very low energies. In Figs. 2. f and 3.f we show the data 
from the Geneva-Saclay group J33| (solid triangles) and 
the very recent, and more precise, data from E865 col- 
laboration at Brookhaven [Q (solid squares). 

Channel (I,J) = (1/2,1) Here the possible states are 
Ktt and K-q. We have plotted in Figs 2.g and 3.g data 



from the following experiments: |35[] (solid square) |3q ] 
(solid triangle). Note that the first set is systematically 
lower than the second, which is newer and more precise. 
Nevertheless, they are compatible thanks mostly to the 
large error bars on the first set. 

Channel (I,J) = (l/2,0) Here the states are also Ktt 
and Ki]. The data in Figs 2.h and 3.h come from the fol- 
lowing experiments: |?5| (solid square), []37f (open trian- 
gle), p|] (open diamond), |3(| (solid triangle), ]39| (open 
squarejT It can be easily noticed that not all the data sets 
are compatible within errors, but once again they can be 
reconciled by assuming a systematic error of the order of 
a few percent. 

Channel (I,J) = (3/2,0) The only state here is irK. 
In this case we have plotted in Figs. 2. i and 3.i, data sets 
from |}6| (solid triangle) and (solid square). The 
latter is somewhat lower than the former, although they 
are compatible mostly due to the large errors in |3q ] 

Channel (I, J) = (1,0) The possible states for this case 
are n-q and KK. We have plotted in Figs 2.j and 3.j the 
ttt] effective mass distribution from the pp — > p(r)TT + TT~)p 
reaction studied by the WA76 Collaboration ffl]. In or- 
der to reproduce this data, we are using 



dE r „ 



cPirri\Ti2\ 2 + background. 



(33) 



Where the c factor accounts for the normalization of the 
mass distribution and the dashed curve in those figures 
corresponds to a background due to other resonances 
apart from the ao(980) (see [[ll| for details). 

Once we have described the data on the different chan- 
nels, we will first compare with the IAM "predictions" 
from the present values of the ChPT low energy con- 
stants, and later we will fit these data by means of the 
IAM. 



VI. THE IAM WITH PRESENT LOW ENERGY 
CONSTANT DETERMINATIONS 

In this section we will comment on the results of ap- 
plying the coupled channel IAM using the low energy 
constants from standard ChPT. Since the values of these 
constants have been determined from low energy data 
or large N c arguments, the high energy results could be 
considered as predictions of the IAM. For our calcula- 
tions we have used: U = 92.4 MeV, M r = 139.57 MeV, 
M K = 495.7 MeV and M n = 547.45 MeV. 

In the second column of Table I we list the values ob- 
tained from a very recent and precise two-loop 0(p e ) 
analysis of Kn decays . Note that the errors are only 
statistical. In the next column we list the central values 
of the same analysis but only at 0(p 4 ). In the fourth col- 
umn we list the values from another set where L\, L2, L3 
are taken from an overall fit to K e 4 and tttt data |43} | 
and the rest are taken from fl. Note that all of them 
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are pretty compatible and, except for L5, the size of the 
error bars is comparable. 





K.i^ decays 0{p^) 


K14 decays 0(p 4 ) 


ChPT 


L\{M„) 


0.53 ±0.25 


0.46 


0.4 ±0.3 


L r 2 (M p ) 


0.71 ±0.27 


1.49 


1.35 ±0.3 


L 3 


-2.72 ± 1.12 


-3.18 


-3.5 ±1.1 


Ll(M p ) 








-0.3 ±0.5 


Ll{M p ) 


0.91 ±0.15 


1.46 


1.4 ±0.5 


Ll(M p ) 








-0.2 ±0.3 


L 7 


-0.32 ±0.15 


-0.49 


-0.4 ±0.2 


Ll{M p ) 


0.62 ±0.2 


1.00 


0.9 ±0.3 



TABLE I. Different sets of chiral parameters xlO 3 . The 
second and third columns come from an 0(p 6 ) and 0(p 4 ) 
analysis of K14 decays [42] , respectively. Note that L\ and L% 
are set to zero. In the third column L\, L%, L3 are taken from 
[43] and the rest from [2] (L4 and Lg are estimated from the 
Zweig rule). 

In Fig. 2 we show the results of the IAM with the values 
given in the fourth column of Table |]. The solid curve cor- 
responds to the central values, whereas the shaded areas 
cover the uncertainty due to the error on the parameters. 
They have been obtained with a Monte-Carlo Gaussian 
sampling of 1000 choices of low energy constants for each 
y/s, assuming the errors are uncorrelated. It is worth 
noticing that these error bands are so wide that the re- 
sults for the other columns in Table || are rather similar, 
even for the central values. Qualitatively all of them look 
the same. 

It can be noticed that the IAM results, even with the 
low energy parameters from standard ChPT, already pro- 
vide distinct resonant shapes of the p, Jo (980), K* , and 
ao(980) (see Figs. 2. a, 2.b, 2.g, and 2.j, respectively). In 
addition, the IAM also provides two other extremely wide 
structures in the (0, 0) tttt and (1/2, 0) nK scattering am- 
plitudes. They correspond to the a (or / (400 - 1200)) 
and k (see Figs. 2.b and 2.h). These structures are too 
wide to be considered as Breit-Wigner resonances, but 
they are responsible for the relatively high values of the 
phase shifts (the strength of the interaction) already near 
threshold. In the last years there has been a considerable 
discussion about the existence and properties of these 
two states (for references, see the scalar meson review in 
the PDG Since ChPT does not deal directly with 

quarks and gluons, it is very difficult to make any con- 
clusive statement about the spectroscopic nature of these 
states (whether they are qq, four quarks states, meson 
molecules, etc..) unless we make additional assumptions 
, which would then spoil much of the model inde- 
pendency of our approach, which is based just on chiral 
symmetry and unitarity. Nevertheless, the simplicity and 
remarkable results of this method gives a strong support, 
from the theoretical side, for the existence of both the a 
and the n. From previous works, it is known that the 
ChPT amplitudes unitarized with the IAM generate the 
poles in the second Riemann sheet associated with the 



cr and the n around ^/s^T ~ 440 - i225 MeV [|j9) and 
yfSp i e ~ 770 — i250MeV |J, respectively. (Let us re- 
member that since these states are very wide, the familiar 
relations M ~ Re ^Js po i e and T ~ — 21m ^/s po i e are very 
crude approximations). We have checked that similar re- 
sults are obtained for the amplitudes of this work. These 
values have to be considered as estimates, since the un- 
certainties must be rather big, taking into account that 
the data in this channels are very conflictive (see Figs. 
2 and 3). The fact that we are able to reproduce these 
states with parameters compatible with previous deter- 
minations is also a strong support for their masses and 
widths, which are in agreement with recent experimental 
determinations both for the a and the k [Q . 

To summarize, we have just shown how the present 
status of both the experimental data and the Li deter- 
minations allows for a use of the IAM despite the approx- 
imations made on its derivation, like the poor description 
of the left cut commented above. 



VII. IAM FIT TO THE SCATTERING DATA 

Once we have seen that the IAM already describes the 
basic features of meson-meson scattering, we can pro- 
ceed to fit the data in order to obtain a more accurate 
description. For that purpose we have used the MINUIT 
function minimization and error analysis routine from the 
CERN program library Q . 

Our results are presented in Fig. 3, whose different 
curves and bands can be understood as follows: As we 
have already seen when commenting the experiments in 
the previous section, and as it can be noticed in Figs. 2 
and 3, there are several incompatible sets of data for 
some channels. In the literature, this is usually solved 
by adding an extra systematic error until these values 
are compatible. We have made three fits by adding a 
1%, 3% and 5% errors to the data on each channel. The 
continuous line corresponds to the 3% case and the re- 
sulting Li values are listed in the second column of Table 
|| . The shaded areas have been obtained again from a 
Monte-Carlo sampling using the Li uncertainties given 
by MINUIT for this fit, which are listed on the third 
column of Table H. Let us remark that there would be 
almost no difference to the naked eye if we showed the 
fit with a 1% or a 5%, neither in the central continuous 
line nor in the shaded bands. Furthermore, the % 2 /d.o./ 
for any of these fits is always 0(1). 

However, although the curves may almost remain un- 
changed when fitting with a different global systematic 
error, the values of the Li come out somewhat different 
from each fit. This is an additional source of error on 
the Li parameters, listed on the fourth column of Table 
II. It can be noticed that it dominates the uncertainty 
on the Li. For illustration, the area between the dotted 
lines in Fig. 3 corresponds to a Gaussian sampling of the 
chiral parameters with the two sources of error added in 







quadrature. 

By comparing the L\ from the IAM fit in Table II with 
those of previous ChPT determinations (in Table I), we 
see that there is a perfect agreement between them. This 
comparison of the complete IAM fit parameters is only 
possible now that we have the full 0(p 4 ) amplitudes, 
given in Appendix B, which are regularized and renor- 
malized following the same scheme as standard ChPT. 
In particular, the agreement in the value of L-j indicates 
that we are including the effects of the rf consistently. 





Fit+crrors 
(curve in Fig. 3) 


MINUIT error 
(band in Fig. 3) 


From data 
systematic error 


L\{M P ) 


0.56 ±0.10 


±0.008 


±0.10 


L r 2 {M p ) 


1.21 ±0.10 


±0.001 


±0.10 


L 3 


-2.79 ±0.14 


±0.02 


±0.12 


Ll{M p ) 


-0.36 ±0.17 


±0.02 


±0.17 


Ll{M p ) 


1.4 ±0.5 


±0.02 


±0.5 


Ll{M p ) 


0.07 ±0.08 


±0.03 


±0.08 


L 7 


-0.44 ±0.15 


±0.003 


±0.15 


Ll(M p ) 


0.78 ±0.18 


±0.02 


±0.18 



TABLE II. Low energy constants obtained from an IAM fit 
to the meson-meson scattering data. The errors listed in the 
second column are obtained by adding in quadrature those of 
columns 3 and 4. 

The threshold parameters (scattering lengths and slope 
parameters) obtained with the IAM are given in Table III 
for the low energy constants in the second column in Ta- 
ble P The errors in Table Eld are obtained by a Gaussian 
sampling of the above low energy constants. Note that 
the experimental values of the threshold parameters have 
not been used as input in the fit, and the numbers we give 
are therefore predictions of the IAM. As we have antici- 
pated before and Table [ill] chows clearly, we are able to 
reproduce the low energy behavior with great accuracy. 
Let us then comment, for each different channel, on the 
results of the IAM fit: 

Channel (I,J) = (1,1) The most striking feature of 
this channel is the p(770) resonance, which, as it can 
be seen in Fig. 3a, can be fitted with a great precision. 
This had already been achieved at 0(p 4 ) both with the 
single H and the coupled (ll]] channel formalisms. How- 
ever, this is now achieved in a simultaneous fit with all 
the other channels, but since we are using the complete 
0(p 4 ) expressions we have a good description of the high 
energy data without spoiling the scattering lengths listed 



in Table III 



This channel depends very strongly on 2L\ ± L3 — L\, 
and this combination can thus be fitted with great ac- 
curacy. The mass and width from a clear Breit-Wigner 
resonance can be obtained from the phase shift by means 
of 



Sij(Mr) = 90° 



1 



d6 



1,1 



ds 



(34) 



s=Ml 



T p = 135.5igQ MeV , in perfect agreement with the val- 
ues given in the PDG [}45f . The errors correspond to a 
Gaussian sampling with the central values quoted in the 
second column of Table [n| and the MINUIT errors of the 
fit. 

Finally, and just for illustration, the inelasticity pre- 
diction from the IAM is shown in Fig. 3d. Note that the 
data values are so small and the claimed precision is so 
tiny that any other effect not considered in this work (like 
the 4-7T intermediate state) would yield a contribution be- 
yond the precision we can expect to reach with the IAM. 
That is why they have been excluded from the fit. 

Channel (I,J) = (0,0) There are three independent 
observables in this channel with data. Concerning the 
7T7T scattering phase shift, plotted in Fig. 3b, we can re- 
produce two resonant structures. First, there is the a (or 
/o(400 — 1200)), which corresponds to a broad bump in 
the phase shift, that gets as high as 50° not very far from 
threshold. This is not a narrow Breit Wigner resonance. 
Indeed it was shown in the IAM with just one channel [||, 
that it is possible to find an associated pole in the second 
Riemann sheet, quite far from the real axis. Second, we 
can nicely reproduce the shape of the /o(980) which cor- 
responds to a narrow Breit-Wigner resonance although 
over a background phase provided by the cr, so that its 
mass and width cannot be read directly from Eq . (|34|) . 

Once more, it can be seen that the scattering lengths 
can also be reproduced simultaneously with the high en- 
ergy data. 

The next observable is the irir — ► KK phase shift, 
Fig. 3c, which can also be fitted neatly. Since we have in- 
cluded the 7777 intermediate state, the fit is somewhat bet- 
ter than with just two channels above the two rj thresh- 
old, as it was suggested in pj|| , but not as much as ex- 
pected (this could be due to our crude treatment of r\ — 77' 
mixing, that we commented at the end of section II). 

Finally, in Fig. 3d, we show the inelasticity in the (0, 0) 
channel. These are the most controversial sets of data, 
since there is a strong disagreement between several ex- 
periments (up to a factor of 2 in the overall normaliza- 
tion), as we have mentioned when commenting the data 
on this observable. 

Channel (I,J) = (2,0) We have plotted the results in 
Fig.3.e. Since only the 7T7r state can have these quantum 
numbers, we are simply reproducing the single channel 
IAM formalism that already gave a very good descrip- 
tion of this non-resonant channel 0|. Nevertheless, let us 
remark that it is now fitted simultaneously with all the 
other channels, and the value of the scattering length ob- 
tained from our fit is compatible with the experimental 
result and standard ChPT, see Table III. 

In addition, once we have a description of this and 
the (0, 0) channel, we can obtain the phase of the e' pa- 
rameter which measures direct CP violation in K — > tttt 
decays Eel]. It is defined, in degrees, as follows: 



e') = 90° - Ofo, - S 20 ) s 



(35) 



For the (1,1) case we obtain M p = 775.7+3;g MeV and 
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Our result is <fi(e') = 38 ±0.3, where the error is obtained 
from a Gaussian sampling of the parameters listed in col- 
umn 2 of Table II with the MINUIT errors on the third 
column. This is in very good agreement with the exper- 
imentally observed value of 4>(e') — 43.5 ± 7. Standard 
ChPT g§ predicts 45 ± 6. 

Low energy Ku decay data. There is no real im- 
provement in the description of these low energy data in 
Fig.3f compared to ChPT, since standard ChPT is work- 
ing very well at these energies. However, these very pre- 
cise data at so low energies ensure that the parameters 
of our fit cannot be too different from those of standard 
ChPT. In addition, they are extremely important in the 
determination of the scattering lengths, in particular, of 
the controversial aoo- 

Channel (I,J) = (1/2,1) As it happened in the (1, 1) 
channel with the p, this channel is dominated by the 
K*(892). This is a distinct Breit-Wigner resonance that 
can be fitted very accurately with the IAM, see Fig.3g. 
From Eq.(§4|) we find M K , = 889 ± 5 MeV and F K * = 
46 ± 13 MeV, in fairly good agreement with the PDG 
p5[ . The errors have been obtained in the same way as 
for the p resonance in the (1,1) channel. 

Channel (I,J) = (l/2,0) Due to the wide dispersion 
of experimental results, our fit yields a wide error band 
for this channel, as it can be seen in Fig.3h. Nevertheless, 
as it happened in the (0, 0) channel, the phase shift is 
of the order of 50° not far for threshold, due to a wide 
bump similar to the a in that channel. Here such a broad 
structure has been identified by different experimental 
and theoretical analysis ||^|||5l],[l6| as the k although 
there is still a controversy about its existence and origin 
p2fl , as it also happened with the a. It is very similar 
to the cr, and hence it cannot be interpreted as a Brcit- 
Wigner narrow resonance. 

We also give in Table III the value for the scattering 
length of this channel, in good agreement with the exper- 
imental data, which nevertheless is not very well known. 

Channel (I,J) = (3/2,0) Since only ttK can have 
these quantum numbers, this is once more the IAM with 
a single channel, which already provided a very good de- 
scription J8|. We show in Fig.3i the results of the global 
fit for this channel, as well as the corresponding scatter- 
ing length in Table III. 

Channel (I,J) = (1,0) In our global fit, the data in 
this channel, see Fig.3.j, are reproduced using Eq.(33). 



The shape of the clq (980) is neatly reproduced in the mass 
distribution. In order to compare the value of the nor- 
malization constant c with experiment, we also show in 
Fig. 4 the result of applying the IAM with the parameters 
obtained from our fit to the experimental data obtained 
from R'-p -► S+(1385)7r?7 and K~p -> E+ (1385) KK 
p3] |. These data have not been included in our fit since 
they do not have error bars, but it can be seen that the 
IAM provides a good description. Once again we are 
using a formula like Eq.(|33|), but with a constant differ- 
ent from that for Fig.3.j and no background. Our result 
is c = 63 ± 15pb/GeV, to be compared with the values 



quoted in |53] where c was taken from 73 to 165 pb/GeV. 

Channel (I,J) = (0,1) Finally, we show in Fig. 5 the 
results for the modulus of the amplitude in the (0, 1) 
channel. In this case, there is only one meson-meson 
scattering channel, namely KK — > KK. Therefore, we 
can only apply the single channel IAM, and in so doing 
we find a pole at approximately 935 MeV in the real axis. 
The width of this resonance is zero, since within our ap- 
proach it can only couple to KK and its mass is below 
the two kaon threshold. One is tempted to identify this 
resonance with the 0(1020) meson, but in fact it can only 
be related to its octet part ui$. The reason is that the 
singlet part w\ is SU(3) symmetric and it does not cou- 
ple to two mesons since their spatial function has to be 
antisymmetric. Consequently we can only associate the 
resonance obtained with the IAM to the octet uig, . 
The position of the pole seems consistent with an inter- 
mediate mass between the 0(1020) and the w(770). This 
state had also been found when using the IAM with the 
incomplete chiral amplitudes 0, and it had been used 
later to study the — > tttt decay within a chiral unitary 
approach The fact that we find it here again con- 
firms that it is not an artifact of the approximations used 
in P] . In addition although the amplitudes used here are 
complete up to 0(p A ) and the fit is rather different, it 
appears almost at the same place, which supports the 
soundness of the results in M . 





Experiment 


IAM fit 


ChPT 
0(p 4 ) 


ChPT 

o( P 6 ) 


aoo 


0.26 ±0.05 


n 9O1+0.0U3 
u.zo±_ 006 


0.20 


0.219±0.005 


boo 


0.25 ±0.03 


0.30± 0.01 


0.26 


0.279±0.011 


«2 


-0.028±0.012 


-0 0411 + 0009 


-0.042 


-0.042±0.01 


&20 


-0.082±0.008 


-0.074±0.001 


-0.070 


-0.0756±0.0021 


an 


0.038±0.002 


0.0377±0.0007 


0.037 


0.0378±0.0021 


«l/20 


0.13. ..0.24 


n+ 06 


0.17 




a 3/2 


-0.13.. .-0.05 


-0 049+ 002 
u.u^_ 003 


-0.5 




«1/21 


0.017.. .0.018 


0.016±0.002 


0.014 




aio 




o.i5i° ;?I 


0.0072 





TABLE III. Scattering lengths au and slope parameters 
bu for different meson-meson scattering channels. The ex- 
perimental data come from [10,55], the one loop results from 
[5,8,10] and those at two loops from [42]. We are using the 
definitions and conventions given in those references. Let us 
remark that our one-loop IAM results are closer to those of 
two-loop ChPT, although the IAM depends on much less pa- 
rameters than the 0(p 6 ) ChPT. 

Finally, we have also added in Fig. 3 a dashed line that 
corresponds to the result with the central values of the 
parameters in the second column of Table II, but where 
we have used the one-channel IAM at energies where 
there is only one state available, the two-channel IAM 
when there are two, etc... As we commented at the end 
of section [V , this approach ensures exact unitarity at all 



energies, but we can see that it generates a discontinuity 
at each threshold. The results are compatible within the 
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wider error bands with the previous IAM fit (the space 
between dotted lines). This was expected since, as we 
have already commented, the difference between the two 
approaches is of the order of a few percent, which is also 
the order of magnitude of the systematic error added to 
the data for the fit. Of course, it is possible to obtain 
also a fit with this method, as it was done in [g and 
the resulting parameters are still compatible with those 
listed in Table III. 



ticular attention to the uncertainties an errors in our de- 
scription, which have been estimated with Monte-Carlo 
samplings of the fitted chiral parameters within their re- 
sulting error bars. 

Summarizing, we have extended and completed previ- 
ous analysis of these techniques in the meson sector so 
that we believe that our present work could be useful for 
further phenomenological applications. 



VIII. CONCLUSIONS 

In this work we have completed the calculation of the 
lightest octet meson-meson scattering amplitudes within 
Chiral Perturbation Theory (ChPT) at one loop. We 
have calculated three new amplitudes, rjrj — > rjT), Kr\ — > 
K-q and Kr\ — > Kir but we have also recalculated the 
other five independent amplitudes, checking and revising 
previous results. The full expressions are given in Ap- 
pendix B in a unified notation, using dimensional regular- 
ization and the MS — 1 renormalization scheme, which is 
the usual one within ChPT. All the meson-meson scatter- 
ing partial waves below 1200 MeV, with definite isospin 
/ and angular momentum J, can be expressed in terms 
of those eight amplitudes. 

Since ChPT is a low energy theory, the one loop am- 
plitudes have to be unitarized in order to reach energies 
as high as 1200 MeV (and in particular the two kaon 
threshold) . For that purpose we have applied the coupled 
channel Inverse Amplitude Method, which ensures uni- 
tarity for coupled channels and it is also able to generate 
resonances and their associated poles, without introduc- 
ing any additional parameter. In addition, it respects the 
chiral expansion at low energies, in our case up to 0{p A ). 
Thus, we have shown how it is possible to describe simul- 
taneously the data on the (I, J) = (0,0), (1,1), (2,0), 
(1,0), (1/2,0), (1/2,1), (3/2,0) meson-meson channels 
below 1200 MeV, which correspond to 20 different reac- 
tions. We also describe seven resonant shapes, namely, 
the a, ,o(770), K*(892), k, /o(980), a (980) and the octet 

This description is achieved with values for the low en- 
ergy constants which are perfectly compatible with pre- 
vious determinations obtained using standard ChPT and 
low energy data. This comparison is only possible as far 
as we now have the complete 0(p ) expression for all the 
amplitudes in the standard ChPT scheme. Indeed, with 
the present determinations of standard ChPT, we already 
find the resonance shapes and we obtain the most distinct 
features of each channel, although with big uncertainties 
due to the present knowledge of the chiral parameters. 

Nevertheless, we have performed a fit of our unita- 
rized amplitudes to the meson-meson data and we have 
obtained a very accurate description not only of the reso- 
nance region, but also of the low energy data, and in par- 
ticular of the scattering lengths. We have also paid par- 
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APPENDIX A: USEFUL FORMULAE 

Here we will give the main results and definitions of the 
different functions coming from the one-loop ChPT cal- 
culation. We are following the notation and conventions 
of g. 

When calculating the ChPT amplitudes, the typical 
loop integrals that appear are, on the one hand, the tad- 
pole integral, i.e, the Feynman boson propagator evalu- 
ated at x — 0: 



d d q 



{2n) d q 2 - M 2 



2M?\ 



Ml , M 2 
s_ log — - 

167T 2 jl 2 



(Al) 



where \i is the renormalization scale, i = it, K, rj, and we 
have extracted its divergent part for d — > 4, with A given 
in (||). On the other hand, the integrals coming from 
diagrams (a),(b) and (c) in Figure |l| is: 



1 



[q 2 -M 2 p ][(q-p)Z 



Mr 



(A2) 



where P, Q = tt, K , r\ and whose divergent contribution 
in dimensional regularization can be separated as 



Jpq(s) = Jpq(0) + Jpq(s) + G(d- 4) 



(A3) 



where 
Jpq(0) 

>00 



= -2A- 



1 



32tt 2 



1 1 
16?r 2 A 

A 

s 



m p 

M 2 log^ 
p 



M o 



log 



Ml 



log 



(s - u(s)) 2 - A 2 



(A4) 



12 



and 



A = Mr 



M% 



f M% 
{M P - 



■ M, 



Q) 



s - (M P - Mi 



Q) 



For the case of a single mass Mp = Mq, the above 
integrals read 



Jpp(s) 



Jpp(s) = 



-2A 



16tt 2 



1 



16tt 2 



1 



, m p\ , 
log— f ) +J PP (s) 



a(s) lo; 



a(s) - 1 



with 



cr( s ) = (1 - AMp/s) 



.1/2 



(A5) 



(A6) 



Note that the above integrals have the correct unitarity 
structure in the right cut which extends on the real axis 
from s = (Mp+Mq) 2 to infinity. In fact, all the integrals 
appearing to one loop in ChPT can be expressed in terms 
of the tadpole and J integrals above fl. However, it is 
customary to express the results also in terms of: 



Jp Q {s) = J{s)-sJ'(0) 



where, from (A4) one has 
1 



J'(0) 



32tt 2 



S M p Mq 
A 2 + 



A :i 



, M Q 



(A7) 



(A8) 



From the above definitions it is easy to check that the 
functions J(s)/s and J(s)/s 2 have well-defined limits as 
s 0. 
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FIG. 2. The curves represent the result of applying the coupled channel IAM using the determination of the ChPT low 
energy constant given in the fourth column of Table 1. The shaded area covers the uncertainty due to the errors in those 
determinations (assuming they were totally uncorrelated) . 
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FIG. 3. The curves represent the result of the coupled channel IAM fit to meson-meson scattering observables which is 
described in the text. The shaded area covers only the uncertainty due to the statistical errors in the Li parameters obtained 
from MINUIT (assuming they were uncorrelated) . The area between dotted lines corresponds to the error bands including in 
the Li the systematic error added to the data (see text for details). Finally, the dashed line corresponds to the use of the one 
channel IAM when only one channel is accessible, but keeping the same parameters of the previous fit. 
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FIG. 4. We show the effective mass distributions of the 
two mesons in the final state of K~p — » E + (1385)-7r?7 (top) 
and K~p Y, + (1385)KK (bottom), the data come from 
[53]. The curves and bands are as in Fig. 3 
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FIG. 5. We show the modulus of the (I, J) = (0,1) 
KK — > KK amplitude. The pole around 935 can be iden- 
tified with the octet ujs (see text for details). Although that 
cannot be shown in a plot, the modulus of the amplitude ac- 
tually becomes infinite. 
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APPENDIX B: ONE LOOP AMPLITUDES FROM CHPT 



Here we list the expressions of the eight independent meson-meson scattering amplitudes to one loop in ChPT. We 
have carefully checked the scale independence and perturbative exact unitary (see section ID). Note that we have 
used eq.(||) to write all the fx and f v in terms of /„■, L\ and Lg, in order to ensure "exact" perturbative unitarity, 
eq.(|2^). Let us first give the three amplitudes that had never appeared in the literature in any form: 

For r\r\ — » r\r\: 



T(s,t,u) 



UK 



-2- {81 [t 2 -su- AtM 2 ] + 14M* - A8M 2 M 2 



K 



378Af*} 



., r , vf2 {M$ - 8M 2 M 2 + 2AM*} + j^(2L\ + 2L r 2 + L 3 ) {s 2 +t 2 + u 2 - AM*} 



(Bl) 



— {\2M*L\ + (3M* - 10M 2 M 2 + l3M*)L r 5 - 36M*L r 6 - 2A(M* - 3M 2 M 2 + 2M*)L 7 



6L r s (2M% - 6M 2 M 2 + 7M*)} 



1 



1 



1927T2/ 4 
ri 



{27 (t 2 -su- AtM 2 ) + 16 (23A/4 - 22M 2 K M 2 + 10M&)} 



+ ^4 |^( 16M i - 1MlfJ m {s) + M$J^{s) + i(9a - 2M 2 - &M 2 ) 2 J KK (s) + [s^t] 



For K°r] -> K°r) 

9t-6M 2 -2M 2 2LI 
T{s,t,u) = — 



[3M* + 12M* + M 2 (5M 2 - to)] 



m 2 3/4 

l - {2 (12LJ + 5LJ) (2M& - t) (2A/ 2 - t) + (12L£ + L 3 ) [(s - M 2 K - M 2 ) 2 + (u- M% - M 2 f] } 



(B2) 



{8{L r 6 - L\)M 2 K M 2 + 2L 7 (M* - AM 2 M 2 + 3M%) + L r s (M* - 3M 2 M 2 + 6M%) + 2L r A t{M 2 + M 2 K )} 

Jit 



/<tt 



48/2 (M 2 K - M») 



{2M 2 K [26M 2 + 69i] - 84M% + 3 [16M* - 50tM 2 + (s - uf] } 



12pM 2 K [m 2 k - M 2 ) ^ 92M ^ " 81M ^ " 6 ° M ^ ^ + M ^ + 18M ^ ~ 2su + 6tM v + 8M v)} 



144/2M2 (Mi - Ml 
2 [153 i - 166M|-] + 144M^} 



{lUtM& - 128M 6 K + [270 



486iM| + 428Af| r ] M 2 



{116M 4 + M 2 K [18AM 2 - 153 1] 



- 9 [lOi 2 + 2su - 3tM 2 + AM*] } 



2304/47T 2 

tJ KK {t) (to - 2M 2 - 6M 2 ) J m {t) (to - 2M 2 - 6M 2 ) (l6M| - 7M 2 ) 



16/4 



216/4 



tJ^(t)M 2 



32/4 



Jkti(s) 



[27s(s - u) + 189MI + 8M 4 + 54uA'/ 2 + 45Af* + 12Af 2 (3s - 2M 2 ) 



18M| (6s -3u + AM 2 + 9M 2 )] 



JKn(s) 



[27s (s - «) + 29M| + 11 A/ 4 + 18A/4 



2M| r (18s + 27u - 47Af 2 - 78M 2 ) + 6Af 2 (9m - 6s + 8M 2 )] 



Jktt(s) 



[Mk (3u + UM 2 - 8M 2 ) 



2M% - 2M 2 K Ml (3u + 5M 2 + AM 2 ) + M 2 (6M% + M 2 (3u + AM 2 ))] 



(M 2 K - M 2 f (AMI ~ 18M| - 6M 2 - 3u) 



6(M 2 K -M 2 ) 2 



J K ,{s)(M 2 K -Ml 



J Kr ,(s)(M 2 K -M v 



J2\ 



+ s <-> u\ 



For K°r] -> ^tt : 
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, 8M% + 3Af 2 + Af 2 — 9t n 22 

T(s, t, u) = - ^= + F ^- {27s 2 + 18su + 27m 2 + VJUMl (B3) 

12^3/2 48V3/2(M|.-M2) 1 V 7 

- 292Afi + 12(5M| - 6i)Af 2 - 32M*) = |9t 2 M 2 + 24Af|- 

+ AMUHMl - 15t) + 2M& [9(s - u) 2 + UMl - 22Mf\ } = ^- |3(s - u) 2 

KV ' " " J/ 16V3/ 2 (Af|--M2) 1 V ; 

2 i i n ?l/f2\ / 7i/r2 ofl/f2M i ' f<V<">„ i ...A/n i o„.A ini7i,f4 oo+fi,f2 ii;i#4 



'/ 



+ 2(3t - 14M£ + 10M 2 )(M 2 - 2M 2 )} + 4 {2(2s + «)(» + 2u) - 192M£ - 23tAf 2 - 16M, 

+ 5M| (13< + 24M 2 ) } - -p 3 - {s 2 + 4su + u 2 - 30M| - 2<Af 2 + 2M* + 6M&(* + 2Af 2 )} 

v3/^ 

+ {3M* [L5 - 2(2L 7 + LI)] + M* [6(2L 7 + L r 8 ) - L\\ - 6L r 5 M*(t - Af 2 )} 

9t - 8Ml - Ml - 3Af 2 r _ , - -, 1 
M4Vf^ [3*^(*) + 4M *-W*)] + 288Vl/4 {-W*) [27s(^ - s) - 45M* + 14M T 4 

- 6Af 2 (9m + 7Af 2 ) - 9Af^ + Aff (36s - 54u + 22M 2 + 156M 2 )] + 3 J K ^(s) [29Af| + 7Af^ 
+ 3s (9s + 3u - 4Af 2 ) - 2A/| (l6s + 9m - 18Af 2 + 3Af 2 ) - Af 2 (40s + 18u - 30M 2 )] 

+ 9 Jk ^ (M 2 k - A/ 2 ) [lOAf£ + 2M$ - A/ 2 (3u + 8Af 2 ) + Af| (3u - 12A/ 2 + 8M 2 )] 
+ (Ml - Ml) 2 (in - 2M% + 2Af 2 ) - b AAjM ( M | - Af 2 ) (Af 2 - - A/ 2 ) 3 

- M J ^ (S) (Ml - Mlf (Ml - Af 2 ) + [, ~ u] j . 

Apart from the above three amplitudes, we have recalculated the other independent five. The reason is threefold. 
First, we wanted them to satisfy exact perturbative unitarity to apply the simplest IAM formulae. This was not 
the case of all the calculations in the literature, even when considering the one-channel case. Second, there have 
been several unfortunate misprints and errata in the published formulae (including some errata made by one of the 
authors). Finally we would like to have a self-contained description of the one-loop calculation, together with all 
the resulting formulae. Nevertheless, when compared with previous analysis, our results are not exactly the same 
because we have chosen to express the amplitudes in terms of only one physical decay constant / w , and we have only 
used the Gell-Mann-Okubo relation to simplify masses if it did not affected the exact perturbative unitarity relation. 
Apart from corrections, the differences are 0(p 6 ). The first amplitude to appear in the literature was 7r + 7r~ — > 7r°7r°, 
although in SU(2) However, we have been able to check also with the SU(3) calculation ]Tc| ] . The result, following 
the notation in Appendix A, is: 

T(s t u) = S - M * M?r Us 2 - Atu - 4sAf 2 4- 9Af 4 ) m Is 2 -tu + 2sM 2 \ - (v,a\ 

+ ± {(2L\ + L 3 )(s - 2Af2) 2 + U 2 [(t 2Af2) 2 + („ _ 2A^) 2 ]} 

J IT 

o j\/r2 

' - {{2L\ + L r 5 ) S + 2(2V 6 + LI- 2L\ L r 5 )M 2 } 



f 1 

./ 7T 



' {30(Af 2 - s)s + 21tu - 56Af 4 } + i (f!^f) + M^^M + (fi 2 _ M ^ {s) 



576^/4 " J 2/4 [ 4 9 

+ _^ | ( t -4M^)(2., + »-4M;)^(0 + [((( _ u) _ 2m;(( _ 2 „ + m;)] ^ (() + [t „ „, 

The K + tt + — * K + tt + one loop calculation was first given in juj. It was correct up to 0(p 4 ) but when expressed 
in terms of physical constants it did not satisfy exact perturbative unitarity. One of the authors gave an expression 
satisfying that relation ||, but there was also a typographical error in that reference. Our corrected result, expressed 
just in terms of f v is: 

T3/2 (s , t> u) = M * + ^- s + JL {(4L[ + L3) (t 2Ml) (t - 2Af 2 ) + {2U 2 + L 3 ) (u-M 2 K - Af 2 ) 2 (B5) 

z Jtt Jit k 
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2L r 2 (s - M 2 K - Af 2 ) 2 + AL\ [tMl + M 2 K (t - 4M 2 )] - 2L r 5 M* {s + M 2 K - Ml) 
8 (2L r 6 + LI) M*M% + 24/ 2 M 2 ( ^2_ Mi) i 2M K fa + 5" " ^M 2 K ) (2Af 2 - t) 
[26s 2 + 21su + 25m 2 - 3Aff (a + 5m + 16Mf)] Af 2 + (85s + 53m - 78Af|) M* - 66M^} 



Mir 



{42M| - Af 2 (5s + 4m - 9Af 2 ) (2Af£ - t) + AM% (l2M 2 - 13s - 8m) 



Ml [lis 2 + 12sm + 7m 2 - 3Af 2 (s - u + 8Af 2 )] } + ?2/2M2 7^ _ M a) { 41M ' - 18 (* + u) Af 4 
[36 (s - m) m + 9 (s + 5u) Ml - 59Af 4 ] Af 2 - 3 [9 (5s + m) - 43M 2 ] Af 4 + 81M^} 

' {3 (s - lOt) i - 6su - 3u 2 - 27M|- + M% (30s - 3i + 21m - 34M I) - Mil (t - 2s - 3m + Af 2 ) } 



1152/ 4 7T 2 
J K it{u) 



4 



{ (s - 5m) u + 5M% - 2 (s - 2m) Af 2 + 5Af 4 - 2M| (s - 2u + 5Af 2 ) } 



16/, 

+ ^ {* (« - 2s + 2A4) + (8s + 3t - 8M|) Af 2 - - [2Af 2 + 6Af r) 2 - 9i] 

+ {Mk fa + 3t - 4Af 2 ) + t (4i - s + Af 2 ) - 4Af 4 } + {2Af | (27s + 18m - 74Af 2 - 51M 2 ) 

+ 29Af£ + 38Af 4 - Ml (36m - 48Af 2 ) - 9 (3(s - u)u - 6sM 2 + M*) } + ^jjfi (M 2 K + Ml - s) 2 

- J ^'jil {ZMfc + 10Af 4 Af 2 + 3sAf 4 + Ml (3s - 4Af 2 + lOAf 2 ) - 2M% (2Af 4 + (3s + 4Af 2 ) Af 2 + 3Af 4 ) } 

fM 2 , /r2 ,2 JK.{u)(s-2Ml-2Ml) a (Af| - Af 2 ) 2 Jg» + (Af| - AT 2 ) 2 J K „(u) 

- [M K - M w ) + (M K - Af J ^ . 

The one-loop TT Q rj — > 7r°?7 amplitude was calculated in JhJ. We give here the result expressed in terms of physical 
quantities: 

_ Af 2 (l3Af 4 + 6t (Ml - Af 2 ) - 9Af 2 M 2 ) ^ (Ag - M 4 Af 2 + 4Af 2 Af 4 ) Mr? 
iS " U) ~3/ 2 9/ 2 (Af2-Af2) + 9/^(Af2-Af2) W 

-^ff2" { 20M * (* ~ 3M v) ~ 25M " + 3 [3 (* 2 + *« + « 2 ) + ~ 9 <] I 

+ -i {2 (ij + L3/6) (t - 2Af 2 )(i - 2Af 2 ) + (L$ + L 3 /3) [(s - Af 2 - Af 2 ) 2 + (u- Ml - Af 2 ) 2 ] } 

+ 1 { [t(M I + M 2 ) - AMlM 2 } L\ + 2{2Ll - L r 5 /3)M lM 2 + AL 7 Ml{Ml - Af 2 ) + 2LgAf 4 } 

J TT 

+ * 4 5 ( 77M ^ + M " ( 154M ' - 72t) - 9 [3 (s 2 + sm + m 2 ) + 8tM 2 - 9Af 4 ] } 

1 f Af 2 - t - 

+ 6/4 1 ~t 3m (<)(16M ^ - 7M - 2) + 4 J ^W( 9t ~ 2M " - 6M ") + J - ( 2 * - M ") 

+ |m^( S ) + ^Jkk(s) (9s - 8M| - Af 2 - 3Af 2 ) 2 + [a ~ «]} . 

Finally, the fsTif scattering amplitudes were calculated in |i~l| ] . They were given in a rather different notation from 
the previous ones. Our result is, for K + K~ — » fsr+fC - : 

T C ,(S, t, M) = 2M K ~ U _ [ 5 ( s 2 + st + 4 2) + 6w 2 _ 13wM ^ _ gM 4 ] (B7) 

6/ 2 M^- 

fj, n j , 2 > lls 2 + 8si+llt 2 + 8MAf 2 r -32Af^ 9 (s 2 + < 2 ) + 24uM \ - 6AMl 1 

+ 2ff \H—2M K ) ^ + i 6 (Mi - Af 5) 
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, N f fi4M2 9M2 2 _ 81 (s 2 + t 2 ) - 36 (s + t) M 2 + SMj 9 (s 2 + t 2 ) + 24mM| - 64Mj \ 
+ W ,< [ MM K ~2M 7r -27u ^ + 2(M 2 -M 2 ) j 

+ -i {2L£(u - 2M|) 2 + (2LJ + L 2 ' + L 3 ) [(a - 2Mi) 2 + (t - 2M 2 -) 2 ] - 4LJuM£ - 2^(« - 2M|)M 2 

J 7T 

- 4 ^5 - 2 ( 2L e + L 8)\ m k} H 2304 / 4 tt 2 2/ 1 K ^ Jkk(u) 

+ 28¥F { 6 ° [S (2S + t} + 4UM ^ ~ 8M ^ ] JKK(S) + 2(9S " " M ' " 3M W^3^ 

+ (9s - 2M 2 - 6M 2 ) 2 J m (s) + 3 [s (lis + 4i - 8Af£) - 8 (s + 2t - 4M&) M 2 ] J^s) + [s <-> i] J . 



And for K°K° -» if +if 

2M 2 , - u /xjf 



T„ eti (s, t, u) = — ~2p~~ ~ 12f 2 M 2 K ^ " S?i + " W ^ (S + 16U) + 36M ^^ (B8) 
^ j lis 2 + 4t 2 + 4s (2j + u) - 8 (a + 2t) M| 24 (g - 2t) M 2 K - 9 (.s 2 - 2t 2 ) + 16M* | 

+ 4f 2 \ 5(U " 2M ^ IMS* + 8(M| - M 2 ) j 

+ 12/ 2< p S ^ + i4A ^ 12M 2 + 2(M 2 -M 2 ) J 

+ A /(4LJ + L 3 )(s - 2M 2 K f + 2L r 2 (u - 2M 2 K f + (2L r 2 + L 3 ) (t - 2M 2 K f\ 

J 7T 

+ ± {4LlsM 2 K - 2M* [4L r 4 + LI- 2(2L r e + Ljj)] - - 2M|)M 2 } 

J 7T 

3 (31s 2 + 4su + 16u 2 ) - 4M 2 - (30s + 57m - 80M 2 -) 
2304/> 2 

+ [« («-«) + 4Mi (2M| - t)] - (9s - 8M| - M 2 - 3M 2 ) 2 + ^^(« " 2M|) 2 

+ ^# [« (7, - 4t + 8M|) + 8 (s + 2t - 4M|) M 2 ] + ^ (2s + t - AM 2 K ) (t 4M 2 ) 

+ ^j|(9t - 8M 2 - M 2 - 3M 2 ) 2 + [t( a + 2t) + 4uM 2 K - 8M* ] + ^(9s - 2M 2 - 6M 2 ) 2 . 
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